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ABSTRACT 


Various  aspects  of  reliability  test  programs  when  the  number  of 
items  to  be  tested  equals  or  exceeds  the  number  of  test  chambers 
available  for  testing  and  when  the  underlying  failure  distribution  is 

exponential  are  studied.     Specifically,  two  bivariate  processes  are 

th 
formed  and  characterized;     (i)  the  time  until  the  r       failure  and  the 

accumulated  test  time  of  the  first  r  failures  and  (ii)  the  number  of 
failures  in  the  time  interval  (0,T)  and  the  total  test  time  accumulated 
in  (0,T).     These  are  treated  in  the  two  cases;  replacement  and  non- 
replacement.     Relations  which  exist  between  contractual  specifi- 
cations for  reliability  demonstration  and  the  underlying  distributions 
of  the  test  program  are  discussed.     Maximum  likelihood  estimators 
are  developed  for  the  distributions  when  necessary. 
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i  Number  of  items  available  for  testing 

c  Number  of  test  chambers  available  for  testing 

p  i~c 

T*  The  time  of  the  j"1  failure  measured  from  test  start 

Bj  The  total  test  time  accumulated  up  to  the  j       failure 

B'j'  The  total  test  time  accumulated  up  to  time    y 

Sj  Tj  -  Tj-1 

N^-  Nunnber  of  items  to  have  failed  in  time  t 

f*(s)  Laplace  Transform  of  f(t) 


I.     INTRODUCTION 

A  reliability  program  manager  shall  be  considered  as  one  who 
has  the  task  to  demonstrate,  according  to  contractual  specifications, 
the  reliability  of  a  piece  of  equipment.     Such  a  manager  is  commonly 
faced  with  certain  physical  constraints  such  as  limited  test  facilities, 
number  of  items  available  for  testing,  time  available  to  complete  the 
test  program,  etc  .   The  manager  is  also  faced  with  various  costing 
considerations  such  as  the  cost  of  the  items  to  be  tested,  or  keeping 
the  limited  test  facilities  in  use  for  a  specific  equipment  vice  another, 
etc.     In  view  of  the  preceeding,  the  manager  is  concerned  with  such 
problems  as,  "what  is  the  probability  of  successful  demonstration  of 
the  specified  reliability  throughout  the  test  program,"     "what  is  an 
estimation  of  the  time  required  to  complete  the  test  program,"  and 
"what  is  the  number  of  items  which  will  fail  during  the  testing."     It 
is  the  purpose  of  this  paper  to  provide  mathematical  answers  to 
questions  of  this  type  and  hence  be  an  aid  to  such  a  manager. 

Throughout  this  paper  it  will  be  assumed  that  all  items  to  be 
tested  have  an  exponential  failure  density  (i.e.  f(x)  =  Ae  and   X 

is  the  failure  rate  of  the  item  to  be  tested.)    The  assumption  of 
exponential  failure  density  is  often  made  for  electronic  devices  and 
it  has  been  shown  [Ref .  4]  that  complex  equipments  tend  to  have  an 
exponential  failure  density.     This  assumption  also  leads  to  a  large 
number  of  mathematical  simplifications.     The  extension  to  other 


failure  laws  is  difficult  and  the  present  analysis  can  serve  as  a 
skeleton  for  that  extension.     The  mathematical  developments  in  this 
paper  heavily  rely  upon  the  properties  of  the  Laplace  Transform. 

It  is  a  common  practice  to  write  contractual  specifications  for 
reliability  demonstration  which  call  for  the  demonstration  of  a 
specified  lower  confidence  limit  on  reliability,    Rl ,  with  a  confidence 
coefficient  #  where  the  item  is  required  to  operate  over  a  mission 
time  T.     Lloyd  and  Lipow  [Ref .   4]  show  that  this  specification  gives 
rise  to  the  demonstration  of  an  upper  confidence  limit    \0  for    \  with 
confidence  coefficient     Y  when  the  failure  density  is  exponential. 
This  means  that  the  estimator    A  for     A  ,  the  failure  rate,   must  be 
less  than    A6in  order  to  assert  that  the  specified  reliability  has  been 
de  monstrate  d . 

Truncated  life  testing  [Ref.   4]  is  a  technique  of  experimental 
design  that  economizes  in  terms  of  either  the  number  of  items 
tested  or  the  amount  of  time  spent  testing,  or  both.     This  paper  is 
concerned  with  those  designs  that  terminate  after  either  a 

(1)  fixed  number  of  failures  or 

(2)  fixed  amount  of  time  . 

Other  (sequential)  stopping  rules  are  popular  [Ref.   4],  but  will  not  be 
considered.     The  design  of  the  experiment  will  trade  cost  against  the 
utility  of  success.     The  total  cost  is  the  sum  of  the  number  of  items 
tested  multiplied  by  the  cost  per  item  and  the  amount  of  time  spent 
testing  multiplied  by  the  cost  per  unit  time  (for  use  of  the  test 
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facilities.)    Thus  the  total  cost  is  a  random  number  because  the  time 
spent  testing  for  stopping  rule  (1)  is  a  random  variable  and  the 
number  of  items  tested  is  random  for  stopping  rule  (2).     Thus  the 
total  expected  cost  is  a  function  of    A  and  the  number  of  failures  for 
stopping  rule  (1),  and  is  a  function  of    A  and  the  termination  time  for 
stopping  rule  (2).     Generally  the  expected  cost  will  increase  as  the 
failure  rate     A  decreases.     On  the  other  hand,  the  lower  the  failure 
rate  the  more  desirable  is  the  equipment  and  the  more  likely  it  is 
that  the  contractual  specification  (  \<  \j  can  be  demonstrated.     It 
follows  that  the  joint  distribution  of  (1)  the  estimator  and  the  terminal 
time  and  (2)  the  estimator  and  the  number  of  failures  will  play  a  key 
role  in  balancing  cost  of  testing  against  the  value  of  successful 
completion.     This  latter  value  is  a  utility  type  consideration  and 
must  be  developed  separately  for  each  particular  kind  of  item.     A 
similar  statement  applies  to  the  cost  coefficients.     It  should  be  noted 
that  many  modern  systems  have  high  item  cost  and  high  test  facility 
usage  costs.     Only  the  methodology  for  treating  the  joint  distributions 
mentioned  above  can  be  treated  in  general.     Also,  the  program 
manager  has  the  option  of  terminating  a  test  program  early  if  the 
chances  of  successful  completion  diminish  during  the  course  of 
testing.     Such  a  decision  should  be  based  on  probability  statements 
concerning  the  aforementioned  random  vectors  conditioned  by  the 
current  status  of  affairs.     Thus  an  understanding  of  these  quantities 
is  needed  when  they  are  viewed  as  bivariate  stochastic  processes. 


It  is  assumed  that  the  program  manager  is  constrained  by  a 
limited  number,  c ,  of  test  chambers.     If  the  number  of  items 
available  for  testing,  i,  is  less  than  c,  then  all  items  are  put  on 
test  at  the  same  time  and  testing  takes  place  without  replacement. 
Without  loss  of  generality  this  case  can  be  indexed  by  i  =  c. 
Otherwise  i  >c  and  the  procedure  is  to  test  with  replacement  until 
the  reservoir  of  spares  is  depleted  and  then  continue  to  test  without 
replacement.     The  case  i  =  c  is  treated  in  Chapter  II  and  the  case 
i  >  c  is  treated  in  Chapter  III.     The  joint  distributions  and  the 
marginal  distributions  and  the  moments  for  each  type  stopping  rule 
are  developed.     It  has  been  shown  by  Read  [Ref .   5]  that  the  normal 
approximations  to  these  distributions  are  valid  if  the  number  of  items 
failed  and  the  failure  rate  are  not  too  small.     Explicit  distributions 
will  be  necessary  when  these  quantities  are  small.    An  indication 
for  the  necessary  size  of  these  quantities  for  the  normal 
approximations  to  be  valid  is  also  given  by  Read  [Ref.  5]. 
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II.     THE  NUMBER  OF  ITEMS  TO  BE  TESTED 
DOES  NOT  EXCEED  THE  NUMBER  OF  TEST  CHAMBERS 

A.      STOPPING  RULE:     STOP  AT  r  FAILURES 

In  this  first  case  (i=c  or  p=0),  all  items  are  put  on  test  at  the 
start  of  the  testing  program  and  the  test  is  complete  when  r  of  the 
c  items  have  failed. 

Let  T-j   =  Time  to  first  failure  from  t=0 

T2  =  Time  to  second  failure  from  t=0 


Tr  =  Time  to  r^     failure  from  t=0 


and  S..   =  T1 


S2  -  T2    -T1 


Sp        Tr       Tr_i 
then  S,  has  density  (c-j+1)  X  e  .  "»-  x  >  °-     The  SJ  are 

independent  and  hence  the   Laplace  Transform  of  the  joint  density  of 
the  S.'s  is: 

f*(srt,s2,...,sr)     =      fjJS^ll _  (1) 

j=1       sj  +  (c-j+1)  X 
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It  is  of  general  interest  to  characterize  the  distribution  of  the 
order  statistic.  Let  the  Laplace  Transform  of  the  joint  density  of 
T1 ,  TP, .  .  .  ,  T^      be: 


1»    '  2- 


T      (u,,  «*t>»' 


"■*>        =         £ 


but         T,  =    Z  £ 


,  and  he  nee , 


f*(uM^...ur\  =       E[e. 


--   E  [e 


-  ji 


N  ] 


Therefore,  letting        A\  -  2   Ul<>       yields: 

"J 


J  ll 


(2) 


Inverting  equation  (2)  on  y,.     gives: 


7T      (c'^,u  y 
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and  noting  that  the  Laplace  Transform  of 

€  VCx-c) 

where  lUL'jgijft^   )      is  the  unit  step  function  translated  to     t-|    ,is 

e 


gives  rise  to  the  joint  density  of  T.|  ,  T„,    .  .  .  ,  Tp  as: 


r-i 


-X  Z  ra         -x<c-rM\tv 


for  0   <t-,   <t2      </...      <  tr 

The  marginal  distribution  of  the  time  to  the  j      failure  can 
be  calculated  easily  from  equation  (2). 

It  has  been  shown  [Ref  .2]  that  the  maximum  likelihood  estimator 
for  A   upon  completion  of  the  test  program  is: 

a  r 

X     =-B—  (4) 

r 

where  B     is  the  total  accumulated  test  time  of  all  items, 
r 

(Alternatively,  this  may  be  shown  from  equation  (3)).     Thus  the  joint 
distribution  of    A    and  the  time  to  completion  Tr  is  tantamount  to  the 
joint  distribution  of  Tr  and  Bp.     This  joint  distribution  has  been 
characterized  by  Read  [Ref  .5]  and  for  completeness  the  Laplace 
Transform  and  the  density  is  recorded: 
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where  U(x)  is  the  unit  step  function. 

When  the  progrann  manager  desires  to  mark  progress  at  some 
tinneTduring  the  course  of  testing,  the  data  available  at  that  time 
will  be  a  complete  history  of  the  testing  including  Nf  and  B-x     ,  i.e 
the  current  tally  on  the  number  of  failures  and  the  accumulated  total 
test  time.     Due  to  the  fact  that  the  exponential  distribution  has  no 
memory,  all  data  other  than  N  •y  and  B^  are  irrelevant  in  terms 
of  making  updated  probability  statements  concerning  Tp  and  Bp  . 
It  also  follows  that  the  conditional  density  of  Tr,  Br  given  that 
N~  =nn  (for  m  <  r)  and  B^=b  can  be  characterized  from  equations 
(5)  and  (6)  by  replacing  r  with  (r~m),  c  with  (c~m)  and  by  translating 
Tp  to  T, — m  +T  and  Bp  to  Bp_m  +  b.     Thus: 


f  TV.tiy.    1Ht*YW      ^Sb  ^M"^     =• 


(7) 


-  ""i 


e  71 


J- 1 

and  letting  k  =  r-m  yields: 


U.,+    (Ui-tA^Cc-^-j-i-O 
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■J 

The  univariate  distributions  developed  in  the  same  manner 
from  equation  (5)  are: 


-  >x,  Cc-Y*4-^      r  _>X, 


r-l 


fr,  tou>-  tV)yX  e  L  l-C         J     f  x,>o      (9) 


/  c~^\  k  A  e.  [i-e  J     V(>r    (10) 


4B.Mrt.    ><>'^        £         ,».>o  (11) 

(  v  - 1 }  J 


i L. j        xt>b  (12) 
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The  moments  needed  for  mean  value  analysis  and  for  the  normal 
approximations  developed  by  Read  [Ref  .5]  are: 


A   j.,      Cc-j+i) 


v  I 


r 

*M  -  f 


V**  (  Br  |   N<y  =Kv.j12r  aT]    - 


V-  Wi 


(13) 


Ccw  {TV ,Br  I  nt^  T3^«  T}    =    -L     £  -- 

1  '  Til  J  ,,.  (C-m-O) 
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Now  consider  those  problems  posed  for  the  reliability  program 
manager.     The  manager  can  successfully  demonstrate  reliability  if 
he  can  produce  an  estimator  for    X    ,  say    X     ,  such  that    A  -    A© 
where     X6     is  an  upper  confidence  limit  for      A  .    with  confidence 
coefficient    V    .       Now  upon  completion  of  the  test  program,  the 
manager  will  have  successfully  demonstrated  reliability  if: 


X     =    !1 £    Xo  (14) 


Therefore,  it  is  required  that  the  total  test  time  accumulated  during 
the  test  program  Br  must  be  greater  than  r/^e    where  Br  has  density 
given  by  equation  (11).     Hence,  if  the  program  manager  is  willing  to 
use  a  design  estimate,  say    XD    ,  for   \    ,  he  has  available  to  him, 
the  probability  of  successfully  demonstrating  reliability,  an  estimate 
of  the  total  time  for  the  test  program,  and  various  other  pieces  of 
information  prior  to  the  start  of  the  testing. 

Now  given  that  at  some  time,    y    ,   prior  to  completion  of  the  test, 
the  manager  desires  an  estimator  for  A     based  upon  the  data  which  is 
available  at  ^  ,  it  has  been  shown  [Ref .  2],  that  the  maximum 
likelihood  estimator      A~f    is 

a  Nr 

at  =  — ey- 

A 

and  by  simply  upgrading  the  distributions  in  question  with      \y    , 

the  program  manager  has  a  continuing  estimate  of  the  most  important 

characteristics  of  his  test  program. 
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B.     STOPPING  RULE:     STOP  AT  A  FIXED  TIME  T 

The  distribution  of  the  order  statistic  T1  ,  Tg,    ...   T^  given 

N      =  r  has  qeneral  interest  and  can  be  obtained  from  equation  (3)  as 
T 

follows: 

Prob[NT^    ITV.--X.]     fTl)Tl>ll.lT,   (t./tijW.t^ 


Pv-db 


K"} 


Noting: 


•Prob  I  NTc^   |Tv  --t,]    =     e  (    C,tT 

Prob[NT  =  r]      -*    (^    [|   -   e"^]^  e"X 


and  therefore 


■fTi|Tx.»*'i.TT  I  NT«r    C  Cl  !  tu*  '  "  '  ^   " 


y)  Ay  e     *s> 
[\  -  e        J 


0  <  t1  <    t2<     <  tp  <  T 


(15) 
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Under  a  fixed  time  stopping  rule,  with  i  =  c,  it  has  been  shown 

[Ref  .2]  that  the  maximum  likelihood  estimator  for    A  upon  completion 

of  the  test  program  is: 

a  NT 


A  =  -5 —  c16) 

BT 
and  the  joint  distribution  of  N-r-,  B-p  is  needed  for  experimental 
design  and  decision  making.     Recognizing  that  N-r  is  a  discrete 
random  variable  while  B_p  is  a  continuous  random  variable,  their 
joint  density  will  be  referred  to  as: 


To  develop  this  quantity,  the  following  relationship  is  needed: 


BT    =    BNt    +    (  T  -  TNt  )    (C  -  NT  ) 


and  hence: 


p  (nt  =  r    ,  bt  4  b  ) 


p 


(   Tp   ^    T,     Tp+1     >  T,     Br  +  (T-Tp  )    (c-r)  ^  b    ) 
letting  Jp  =  Bp  +  (T   -  Tp)  (c~r),  equation  (17)  can  be  written  as 
Prob      (tr    k   T,     Tp+1    >   T,     Jp   ^   b) 
and  this  probability  is  equal  to: 
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where  f(x^  ,  Xp,  x^)  is  the  joint  density  of  Tp,  T   _|_i  ,  Jr< 
Let  the  Laplace  Transform  of  this  density  be: 


and  noting  the  following  relations: 

v- 
£    =     I      ^J 

V  +  l 


-,.,  -  1 6, 


4  =  i 

IT 


JV  -  (c-^  i  -  2.  cv-j+o  ^ 


and  hence  make  the  following  transformation  of  variables  in  equation 
(1): 

gives  the  Laplace  Transform  of  f(x-|  ,  X2>  x     )  as: 
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Now  inverting  on  u -|   and  denoting  the  result  Iu  .   f  * ,  yields: 


j^r*  \^  tc*w(^ 


e 


(uj  +  ^v'        (  u^  ^  U-^A) 


J-O  4 

Now  by  factoring  g_  from  under  the 

summation  sign,  and  applying  the  binomial  theorem  to  what  is  left, 

yields: 


'r 


• 


Next  invert  on   U2: 


Iu^f*  = 


W 


(u^Vj 


r-  l 


e-v\cx'lv^  -.f'tfU.-x/i 


.  U,  lt-»\T       -  Xt.U—'I    > 
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Hence 


J  ~oo 


0  T 


I  u-i  u-i.  T       dxu   <i 


-U4'(C-^T        -T(c-r\X 


L    I    -   C  J  (19) 


(  Ks  *  X^ 

which  can  be  rewritten  as: 


r"J  /     \     ^  "  TU-^Uy^ 


£k*  +  X} 


(20) 


Therefore,  upon  performing  the  final  inversion: 
fi 


NT'BT 


(r,b)=     e  X      (^\ 


r-\ 


£     l-lH(r)      [b-TU-^]        Ulk-TCc-j^    (81) 


bid 
which  is  the  desired  joint  density. 

Now  by  considering  the  fact  that 


^>T-}    c    ^ 


-XT  n  r* 


li-c-ATJ 


and  dividing  this  probability  into  equation  (19),  yields: 


E    e 


-u,Br 


JMr--*")    = 


/  -   d. 


-T(u4+\^  n 


Ar  C 


-    UvCC-v^  \ 


L         ^3 


Ua    +  X 
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and   noting 


*I  (*K«-     e     _c L£  -iJ  ,   (22> 


which  can  be  rewritten  as: 


(u3a)  J 

The  first  and  second  moments  are  obtained  from  equation  (22) 


as: 


\Nr]   --    c    [i-  e-*T] 


-AT 


e{Br]    =     C.T    -    i  ( ^  -l+ft"'     J 

.  XT  r  -XT  ") 


(24) 


Vm*  |Br}  =    if  I-  2XTe_A'-  ^"aATl 


->r  r         -xr 


Co,  (Nt.Bt]  -  ^  e"   '  {  '  "  c~      'XTj 
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For  purposes  of  marking  progress  when  the  elapsed  time  of  the 
testing  is   "f<T,  the  memoryless  property  of  the  exponential 
distribution  again  makes  the  updating  of  probability  statements 
relatively  easy.     Given  N  =  m,   B  ^     ='dq,  the  distribution  of 

N-p  and  B      can  be  obtained  by  replacing  r  by  r-m,  c  by  c-m,  T  with 
T   -  X     and  translating  N-p_  v     =  r_m>     B-j-_~   =  x  -  bg. 

At  this  point  the  reliability  program  manager  is  again  armed  with 
all  the  instrumentation  required  to  have  continuing  estimates  of  the 
most  important  characteristics  of  the  test  program.     It  should  be 
noted  for  this  case  however,  that  the  estimator  for    /\    ,  say    )\    ,  is 
now  the  ratio  of  two  random  variables  and  hence  it  is  this  ratio  which 
must  be  less  than  the     A  o  designated  in  the  contractual  specifications, 


i.e.: 
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III.     THE  NUMBER  OF  ITEMS  TO  BE  TESTED 
EXCEEDS  THE  NUMBER  OF  TEST  CHAMBERS  AVAILABLE 

A.     STOPPING  RULE:     STOP  AT  r  FAILURES 

The  development  in  this  chapter  will  generally  follow  the 
development  of  Chapter  II.     Testing  takes  place  with  replacement  as 
long  as  possible  and  testing  is  complete  at  the  r^    failure.    At  the 
start  of  the  testing,  there  are  i  items  available  for  testing  but  only 
c  test  chambers  available,  where  i    >  c .     This  chapter  is  restricted 
to  the  case  where  r  >    p  =  i~c.     To  start,  let: 

T^  =  Time  to  first  failure  from  t  =  0 

To  =  Time  to  second  failure  from  t  =  0 


Tr  =  Time  to  r      failure  from  t  =  0 


and  S-|  =  T1 


S2  =  T9  -  T. 


Sr  -  Tr   -  Tr-1 

-cAx 
then  S:  has  density  c  A    e ,       X>  0  for  j  =  1  ,   2,    .  . .  ,   p  and 


\       -(c+p-j+1 )    A    x 
A   e 


has  density  (c+p-j+1)    A   e  ,  x  >    0  for  j  =p+1  , .  .  .  ,  r 

Hence  the  Laplace  Transform  of  the  joint  density  of  the  S-'s  is: 
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f*(s1  ,   s2,    ...  ,sr)  =  -         ,       I (25) 

Therefore  the  Laplace  Transform  of  the  joint  density  of 

"T"i  »  Tp,  .  .  .  ,  T     becomes: 

TT__^A__     tt         (ctp-j-n  )  ^ 
f*(Ul,   u2,...,up)  =   J  I ^ |(     ~ (26) 


Now  inverting  equation  (26)  on  u-j  yields: 


P         p  1  =  3 


(cA)r  & £ TT  — W 

1-3  K-j 

-pr  (C  +  P-J  +  I  W 


Hence  the  inversion  on  u..  ,  u    ,  .  .  .  ,u     produces: 


J>p*l     (  C+P~  J  +i  )A   4  £    u 


and  0    <:   ti  <,     t0  ....      <  t 

'  2  p 

and  the  format  has  now  become  that  of  equation  (2)  and 


hence  the  joint  density  of  T1  ,  T    ,    .  .  .  ,  T     is 


r-i 
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The  maximum  likelihood  estimator  for     A     Using  a  stopping  rule 
based  on  r  failures  and  assuming  r    >    p  is  developed  in  Appendix  B 


as: 


a  =     r 


Br 

Thus  again  the  joint  distribution  of    A    and  the  time  to  completion  is 
tantamount  to  the  joint  distribution  of  Tr,   B    . 

The  development  of  this  chapter  differs  slightly  from  that  of 
Chapter  II  in  that  the  Laplace  Transform  of  the  joint  density  of  Tr, 
Tm,  Trn+^  ,  B    ,  and  J       (to  be  defined  later)  is  developed  and  from 
this  Laplace  Transform  the  desired  combinations  of  distributions 
are  extracted.     This  method  allows  for  a  clear  presentation  of  the 
relations  existing  between  the  various  random  variables  involved  in 
the  process. 

Consider  the  joint  distribution  of  Tp,   Bp,   I\L,   B  ^ 
and  again  note: 
Prob      (  Tr  ^  t,   Br  ^  b,   B„  -  z,  N  _  =m  )      = 

Prob     (    Tr  6     t,  Tm  £  t       ,  Tm+1  >  X      ,   Bp  ^    b,  Jm   ^    z  ) 

where  J_  =  B_  +  (  X     -  T     )  (i-m)  for    p  <£     m   £    i . 
m  m  m  l 

(The  case  m   ^    p  is  straightforward  and  will  not  be  treated 
explicitely .) 

This  probability  is  equal  to: 

\         dx,      \       dxj.    \       clx3     \      d  xH      \        c\%5     f  (  *i  )*%>  *3  >  *<,  )*6  ) 
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where  f(x*  ....  ,x,-)  is  the  joint  density  of  T    ,  T„,  T      .  .    B^   J 

v    i  '  '    5y  J  J  r'       m»       m+1  ,     r>    m 

Now  let  the  Laplace  Transform  of  this  density  be: 
f*(Ul,...,u5)  = 


and  note  the  following  relations: 


r 


* 
j 


r     Fi     J 


m+1 
T — _i.-\   —  t  S 

j 


m+         i=1      -j 


B 
r 


h  r 

=  c"5ZSi    +     2Z    (c+P"J+1)S. 
j=1     J  j=p+1  J 


P.   ^  m 

Jm  -(i-m)  -T  =  (m-p 


X-  SY  +       T~      (m"J+-1)     si 
j=1  j=p+l  J 


Hence  the  transformation  of  variables 

s.  =  u1  +u2  +   u3  +  cu4  +  (m-p)u5  j=1  ,  .  .  .  ,p 

Sj   =  u1  +  u2  +  u3  +  (c-p-j+1)u4  +  (m-j+1)u  j=p+1  ,  .  .  .  ,m 

s;  =  u.  +  Ug  +  (c+p-m)u4  j=m+1 

Sj  =  u1  +  (c+p-j+1)u4  j=m+2,...,r 
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can  be  applied  to  equation  (28)  and  (noting  that  c+p   =  i),  the  Laplace 
Transform  of  equation  (28)  becomes: 
f*(u1  ,  .  .  .  ,u5)  = 


TT 


U-j+n  A 


jtp-fi         u,-t  u2-+  u3  +  (  u^-f  A)  ( «!-j  +  0  +  U5  (^-j-fi) 


IT 


U,  -f  u3  4  (0^4^  )    (l-/>i  ) 


J-^rl.  U,    4     (  U^4cO     V    L    -T   +    I    ) 

Now  setting  u      =  u     =  u     =    u     -      0      yields  the  Laplace  Transform 


of  density  of  T     to  be: 


TT  (30) 


f*(u,,.;o/0,0)=(_£iL)P+'-]T  U-J.4'.)- 


Equation  (30)  shows  that  T     has  an  Erlang  distribution  with  r 
stages  and  the  indicated  rates.     Now  equation  (30)  can  be  written  as: 
f*  (u<,  0,  0,  0,  0)  = 


(31) 


J       /  i  \  P4i 


([-r)i-         fa        (u.t^V   '  (u,*(t-i+'i)AJ 


Where  72  '    - 


J  tj-f  -2.)]     (r-j)l 
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and  using  theorem  1  of  Appendix  A,  the  density  of  Tr  becomes: 


J*P+«- 


4- 


p  +  I 


I   (j-p-0 

(32) 


H         Z(-() 


(P-j+i)pfi   *- ;  (mo  i 


Similarly,  the  setting  of  u^  ^Jurj  =  uo  =   u5  =  0  in  equation  (29) 

yields  the  Laplace  Transform  of  the  density  of  Br    to  be: 
-*  /         \  \    r 


^0,0,0,^0)     =  (_A-.       ) 


(33) 


which  is  the  Laplace  Transform  of  a  special  Erlang  distribution 

with  r  stages  and  rate   }\     and  hence  the  density  of  EL  is: 

'    ,        vr-i     s-  A  x 

(-  (x)    =     i ,        X    >  o 

O-l)    I  (34) 

Now  considering  equation  (29)  with  U2  =  Ug  =  u,-  =  0,  the 
Laplace  Transform  of  the  joint  density  of  Tp,   Bp  is  produced: 
f*(u1?  0,  0,  u4,  0)  = 


cA  \p+'       -A-  I^J-mW 


0|t^U^  +  A)c     /  J=p  +  2_         U,  +  (  L-j+l)    (    U^  +  /|) 


(35) 
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which  can  be  rewritten  as: 


sT    ~  ^  +  ' 


(/.-^!  Lu.,    *  (*v*X\c"]P*'        (UntX^-p-i 


»-l  y..;-, 


HTj 


5         l-n 


0'.         [  cx.+    (,u4  +^CX,-j'-0  J 


therefore: 

Iui  f*  =  (36) 


Xr  cp+l  U-fr-Oi  )  V         l-J^"' 


5 


U-^>        (uv  +  xr-P-1    ^,    (j-p-OUv-j-.V 


^ ^r     + 


U^^'Cp^*'     f-  (  J     Cp-k\» 


K- O 
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Now  by  inverting  on  u4      the  joint  density  of  Tr  Bp  is  obtained  as: 
Iu1   u4    f*=fTr   ■   Br  Cx1  X4)=  (3?) 


r-i 


p+i  (L-f-i-)i  y    (-i)r"J" 


(r-2)l      (j-P)p  +  ' 


P 


Nr-P+k-2  , 


C-D    e 
K=o  (r-f+k-t)!     (f-l<)!    (f-j)^ 


By  considering  the  logarithm  of  the  transform,  Read  [Ref  .5] 


has  shown: 

E(Tr): 

+ 

i 
A 

r-p-l 

I 

i 

(c-j.) 

Vflk  1  T  1 

( 

P 

1 
t       

< 

1 

1     MIS    y       1   r    f 

j.-o          v 

-i  J1 

(38) 


r 


E  I6r)  =     -j 


(.  o-  v   L    \  v  )  <^  *  /     ~~ 


Pc  I  > 


-r 


A 


XL     j        (c-t) 
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The  program  manager  is  now  armed  with  the  required 
distribution  theory  to  determine  if  the  reliability  testing  has 
successfully  demonstrated  the  specified  contractual  reliability.     If 
the  manager  desires  to  mark  progress  prior  to  the  completion  of  the 
testing,  the  data  available  at  some  time,     I     ,  will  be  a  complete 
history  of  the  testing  including  N<y      ,   B  <y     .     Again  the  memory- 
less  property  of  the  exponential  distribution  makes  N  <~j-    ,   B'y     the 
only  relevant  data  in  terms  of  making  probability  statements 
concerning  the  density  of  Tp,   Br .     Appendix  B  develops  the  maximum 
likelihood  estimator  for    X    ,  say     A^  >  for  a  fixed  time   T    ,  as: 


It  also  follows  that  the  conditional  density  of  Tr,   Br  given 
N~=  m  (for  m  <     r)  and  B  y   =  b  can  be  characterized  from  equations 
(32)  and  (34)  by  replacing  r  with  (r-m),  i  with  (i-m)  and  by 
translating  Tp  to  Tr_m  +  T    ,  and  Br  to     B  +  b,   keeping  track  of 

whether  m  £    p    or    m    >     p.     Thus: 

with  the  appropriate     translations  on  the  indices  depending  on  the 
value  of  m. 
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The  conditional  moments  are  given  by: 


v-p-i 


e{t>  /*T.M«rt}  =  \  *  ?-f?  +   1 .  ?.  ■TvT) 


J 


r-p-i 


VaR?  Ty  I  Nr  ***  Up^*\    -     -^— ■       +     -L      C         — ; (39) 


,T- 


Var^Tv-I  NT«v*C>t^   s       jt    21        Ic^"^ 


J-^-p 


E  f  B^   I    Nr~r»l    BT:lo~j     -     t>    +■ 


V--   Wn 


J  x1- 


Again  the  manager  is  confronted  with  a  maximal  failure  rate 
as  derived  from  the  contractual  specifications  on  reliability.     The 
system  or  equipment  under  test  in  this  case  must  develop  a  total  busy 
time  Br  which  is  greater  than  r/ \0    .     At  this  point  the  manager  is 
equipped  with  the  necessary  managerial  aids  for  decision  making 
throughout  the  test  program. 
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B.     STOPPING  RULE:     STOP  AT  A  FIXED  TIME  T 

Under  a  fixed  time  stopping  rule,  with  i    >   c,  the  maximum 
likelihood  estimator  for     A     upon  completion  of  the  test  program  is: 


\.-5 


and  again  the  joint  distribution  of  Nj,  B-tj- js  required  for  experimental 
design  and  decision  making.  Using  the  same  convention  as  in  Chapter 
II  yields  the  joint  density  on  N-p     B 

9  Hrsb  (  AV=  *«  >   bt  -  \> ) 


'^r-&T 


(«fyb)    - 


3    b 


in  this  case  however  the  procedure  is  not  quite  so  straightforward 
as  it  was  in  Chapter  II.     Equation  (29)  requires  some  understanding 
of  the  indices.     Proper  manipulation  of  the  indices  yields: 


E     e. 


-  U2T,v,  -  U3  TMt|  -  Us    J, 


M 


-  U5-  C  T 


^A 


KVH-  I 


Uj.  4  U3  +  c  A 


A 


u3  *    c^ 


whe  re  m    ^    p . 


It  should  be  noted  that  B-j-  is  a  degenerate  random  variable  for  this 
-Dz7*  ~U3  T„it  -U6   TM 


case . 


E\e 


-Ur     CC-OT 


c\ 


L  ^2.+  u2  ^  \j5  +  c^ 


P+i 


ft-i  U 
u3  -v-  f  L-i)h 


Where  m  =  p  +  1  . 
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£    e. 


-  u-v.    •  *v»  -^Thi^i    -W.J-   J 


*-*\ 


(40) 


-  Utf  (^-vv»^T 


cX 


-\ 


P4I 


Uj,   -t   U^      f     Cwx-p^Ccj-    4    cX 


w\ 


TT 


where  p+1    4     m    <.     i 
The  following  relationships  should  be  noted 


£    e 


fe"u        )nt-v*}« 


Equation  (40)    with  Qg  U3=  0 


-Kr  Sr 


r-o 


-Ui-  Br 


'''}  '    Z   rW  (nt--^  t 


_  f    -u.Nr- ur  Br")  <r     ~  r  ->      -uv"   - 

He  -    £    p<^  {  Mr*"]  C        E 
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The  marginal  distribution  of  N-j-  is  obtainable  by  considering  for 
each    r,r  =  0,    1,2,    ...i      ,     the  probability  N-j-   =  r.       For 
r  =  0,    1,2,     ...p,  the  process  is  a  Poisson  process  with  rate 
Hence: 


-,  ,  v*         -    CAT  (41) 


yl 


r   =   0,1 , . . .   p 
For    r  =  p  +  1 


and  from  equation  (40),  the  Laplace  Transform  of  this  joint  density  is: 


■f      t-u^u^    = 


C  X  I  C  -  ^    A 


Inversion  in  u~  yields: 

Inversion  on  U3  yields: 

P! 
X3  >   XL 
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Now: 


J  dx3  I  u%  Ma  -f*  - 


cp"x?"  xJV^e^"""* 


Pi 


O  <    *,   <    I 


Finally, 


c  -.  P*i     -T(c-iU 


UT) 


(42) 


Now  for    p      p  +  1  ,  the  Laplace  Transform  of  the  joint  density  of 
Tr,  T      ,    .  is  given  by  equation  (40)  as: 


-f^Cu^u.^  - 


[     c  x 

Ux    *      U.^ 

V     ^X    . 

_p+\ 


Y~ 


JSP 


t2- 


u 


U  -  \  +  ^N  a 


(43) 
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Rewriting  gives: 


PH       r+|       (  L-p-l)  ( 

c       A 


i.±ui  — i y  (.or"j 

3^  Pfi 


(L-r-i)l        UyfU-r)*      /_  (  jr-p-i)  I  (r--f)l 


K+U3+c;\   V    '       V  V  u-j+OJ 


I  U-     f  =        C  ?i ; >       (-  I 

-  x2  [u3  +  (L-J-f/)A] 


(i-f-M!  (r-j)!      (ifr.-j)'*'  yJtl 


krtf 


The  only  terms  involving  u     are 


-XZ\J3 

e 


u3+  tl-r)^ 


and  the  inversion  of  this  term  yields 

x1(i.-r)A         -*.»,  C  L-r)A 


e 


e  u  ( x^-x^) 


39 


The  only  ter.m  involving  x3  is  in  the  previous  inversion,  hence 


r 


CO 


J 


In,  w.,   *         d 


iv  ^3 


e  c      A     ( *- 


U  -^ 


J "-  P^ 


2V.  (r-jV. 


->-**.  U-J4^A 

Cj  -  p  -n       A 


I 


P-h  -J  )        A  K" 


P-K 


CP-*V. 


The  first  term  involving  x2  is 


and 


8  dy^    - 


/  -  e 


-T  Ir-jtO  A 


Cv-  \  *»>  A 
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The  next  tenm  is  to  note  for  K   =  p 


*L(r-P)A  I  ,.eT(r"fM   ] 


p e_ 


and  for  all  other  terms  inside  the  sunn 


P-i        ri  ,       J<  p-k        -  **.    (r-p)A 

<c       f  (-i )  xi         £ 

£-        J         (f>  +  l)Ktl     ^kti  cp-  k)    I 


P-I  - 


+  1        \P  +  Z        ,  .     a  _|<  +1 


~        + 


<=D 


(f4    l-J)K*l       ^P.Z        (f_p)    p 


-  t  Cr-p)  ^ 

^+'  £ 


T  r 


l-o 


Therefore  the  probability  that  Nj  =  r  is: 


r  0  U  J 
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-T(i-r^      cPcl   Af 
0--O   1 


i  «■* 


r-J 


i-  p+i 


(j-p-Z)l        (r-j) 


-  e 


-T  (r-j-n)  A 


<re-')p+'AftlM+o   +    (f^-j)ptl   Af+i  (r-ri 


+       (44) 


K+l  -T(r-fU   £k  Tf-*-(- 


l  Ur-ff-^lo*^1  Af+2       (f^-il'V*1    ,Z- 


.(r-f^pH-jr1  A 


(f-k-1)!  [(r-pW] 


m 


ho  ■  -  t=» 

From  the  equations  thus  far  derives  the  moments  of  the  bivariate 

distribution  can  be  developed  as 


Ef*rl- 


E{6T]=- 


I    r   P{VT=r] 

a  u5 

U5-  o 

e  b1 


if  - 


^    a 


(45) 


■r  o 


Now,   in  view  of  the  convention  for  fN  p      (m,b),  this 

bivariate  density  can  be  derived  from  equation  (40),  in  that 


r T  C  °° 

f»T,6T(^b)=     Jo      dX*      JT       ''**     Iu2U3U5      f 


(46) 
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Where  f*  is  the  appropriate  Laplace  Transform  depending  on  the 
value  of  m.     To  start  let  m  £  p  and  hence 


^■%^ri^^     = 


-  usC i     r 


Wt-H 


cA 


c\ 


Ui+    uj  +  tXJ  u^  +  c  X 


Where  e   U5C  '    is  the  Laplace  Transform  of  the  degenerate 
random  variable  B-r.     Therefore: 


vv\  ! 


When  m   =  0,    1  ,   2,    ...   p;  and  b=cT 


Now  for  m  =  p+1 


(47) 


r(^uJ,o  = 


-u^((.-^Tr 


C  \ U 


->U 


U-^X 


p« 


h  >  X*. 
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T 


(u4-  +  /0p+l  ,fr    Lp-^I  (^m<+« 


;»o 


Inverting  on  U5  yields: 
fNT,   BT  CP+1  »  b)  = 

e"^  UX^^    [b-TU-»Q?     Ulb-TCc-,^  (48) 

P! 


It  should  be  noted  that  another  term  would  have  appeared  in  this 
inversion  except  for  limitation  (c-1  )T   £    b  £    cT. 

By  a  development  similar  to  that  of  equation  (48),  it  can  be 
shown: 
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1nt>*t    (**,  b^"= 


e        c   A       c : 


U-v^v, 


r,  _    V^  -  1  r  -iVv--\  (49) 

(  J  -P-0^1    (v*-r^  tvw-0! 


,-i 


Where  m   =p  +  2,   p  +  3,    ...,   i  and  all  terms  in  square  brackets 
are  followed  by  a  unit  step  function  with  the  same  argument. 

The  reliability  program  manager  nowhas  all  the  distribution 
theory  necessary  in  this  case  to  assist  him  in  making  decisions 
about  the  testing  program.     The  formulations  look  somewhat 
unmanageable,  but  with  specific  numbers  these  formulas  should  not 
be  too  difficult  to  solve.     If  the  numbers  become  large,  use  of  the 
normal  approximation  derived  by  Read  [Ref.5]  would  probably  be 
more  practical . 
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APPENDIX  A 


THEOREM  ONE 


Theorem  one  is  stated  as  follows: 


If  f*(s)  is  the  Laplace  Transform  of  f(t)  such  that 


f  *(s)  = 


1 


TT 


(s'+a)  ' '  (s+b) 


»  a  ^  b,  then 


hi 


fCO  = 


-at  £^ 


(a-b)a  O-*)^ 


M 


K*o 


vn  —  i  -  k 
[t  (  b>o()]  (50) 

(n-i  -"k  )  I 


Proof: 


By  the  theorem  of  residues  of  complex  variables, 


fCtU 


s-t 


(s  +  *)n 


cf  '(     S  t  h       ) 


tn.rO  ! 


(51) 


Now  by  the  binomial  expansion  for  derivatives  of  products,     iff 

st 

then  the  (n-1)      derivative  with  respect  to  s  is: 


Jr0 


Now  substituting  equation  (52)  into  equation  (51)  yields: 
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St 


f(t)  = 


(s+a) n 


s=-b 


st  n-1 


i 


(n-1)!  (s+b)n      fi6      C    1KS   DJ  (n-1-j)! 


s=-a 


He  nee : 


-bt                       -at      n-1 
&_    „.n      +       -S— r-    ^    (-1)J(b-a)n-1-Jtn-1-J  1 


'(a-b)n  ^~F    ^ 


(n-1-j)! 
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APPENDIX  B 

MAXIMUM  LIKELIHOOD  ESTIMATORS   FOR  THE 

FAILURE  RATE 

The  estimator  Ay  for    A    upon  completion  of  the  test  program 
(i.e.  upon  obtaining  r  failures)  is  developed  as  follows: 

The  likelihood  function,   L(    X   )»  based  on  the  samples  of  S;,   is: 

P    -cXl*j     tt        t  v  v      -  (C+p-j+OASj 

;  j  *  p + > 

Now  letting  0  (    A    )  =  Ln(    L  (    A   ))>  yields: 

?  * 


J=t 


Jrf*" 


p  V-  v- 

Setting  equation  (53)  equal  to  zero  and  solving  for   X     yields 
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and  recognizing  the  right-hand  side  as  Br  yields  the  estimator  for 
to  be: 

Now  the  maximum  likelihood  estimator  for    A    at  some  fixed 
time    %    ,   prior  to  completion  of  the  testing,  where  N         =m^.p   is 
given  by: 

A 

m 


A  = 


In  the  case  where  N  ^    =  m   >    p,  the  following  data  is 
available:    T    ,  T„,    .  .  .  ,  T      and  B  .     Letting: 

si=Ti 

So  =  To  —  Tm 


<^         =  T"        —  T 
°m        '  m        '  m-1 

and  noting  that  the  probability  of  zero  failures  in  (*C     -Tm),  as  a 
function  of    X    ,  is: 

-(c-(m-p))   A     (x    -T^) 
the  likelihood  function  becomes: 


e  m' 


P         \      -c  A   s      rn 
L-(A)=TTcAe        n       J     TT        (c+p-j+1)    A  . 

j=1  j=p+1 

e-(c+P-j+i)  A  Sj-(C-(,— ))  Act  -Tm) 


49 


Now  letting  0(    \    )    =    Ln(L(  j\     )),  yields: 


0( 


P,  m_ 

>  )    =    pLn(c    }\    )-C>    2,     Sj    +2-.H   Ln((c+p-j+1)  A    )   - 

J-1        J         j=p+.1 


j=p+1 


(c+p-j+1)   X    Sj  -  (c-(m-p))  ^     (t     -Tm) 


P  m  m 

0  ■  c  /\  )  =  -E--  c  y  Sj  +   y    _l_  -  2_    (c+p-j+i  )s 

A  J=1      J  j=p+i     A         j=p+i 


J 


-  (c-(m-p))  (x  -Tm) 


Setting  0'(   ^    )  =  0,  yields: 


h  nn 

HL     =  c     i     S .  =     %  (c+p-j+1  )S      +    (c-(m-p))(X    ~Tm) 

^  3=1       J         j=P+1  J 

Recognizing  the  right-hand  side  as  B/y    yields  the  estimator 
for    n  to  be: 


>^     = 


t 


M 


(55) 
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